Static and Dynamical Phyllotaxis in a Magnetic Cactus. 
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While the statics of many simple physical systems reproduce the striking number-theoretical 
patterns found in the phyllotaxis of living beings, their dynamics reveal unusual excitations: multiple 
classical rotons and a large family of interconverting topological solitons. As we introduce those, 
we also demonstrate experimentally for the first time Levitov's celebrated model for phyllotaxis. 
Applications at different scales and in different areas of physics are proposed and discussed. 
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Phyllotaxis, the study of mathematical regularities in 
plants, challenged Kepler and da Vinci |T], inspired the 
Bravais lattice of crystallography [5] , and may well have 
motivated humanity's first mathematical inquiries [Tj . In 
groundbreaking work, Levitov proposed [3] that the ap- 
pearance of the Fibonacci sequence and golden mean 
in the disposition of spines on a cactus is replicated in 
physics, in the statics of cylindrically constrained, repul- 
sive objects [3]. Here we prove experimentally Levitov's 
model, and describe for the first time the intriguing col- 
lective excitations of the phyllotactic geometry: multi- 
ple classical rotons and a huge family of interconverting 
topological solitons. A simple geometrical mismatch un- 
derpins all these phenomena: nearest neighbors in one 
dimension (i.e. along the axis) are not nearest neighbors 
in the full three dimensions. Being purely geometrical in 
origin, dynamical phyllotaxis could occur in many phys- 
ical systems, including trapped ions in cylindrical po- 
tentials H], Wigner crystals in curved nanostructures, or 
crystalized ion beams [SJ|3]. 

First we review static phyllotaxis as in Levitov [3]: as- 
sume that a set of particles with long range repulsive in- 
teractions, when confined to a cylindrical shell of radius 
R, forms a helix with a fixed angular offset f2 between 
particles and a uniform axial spacing a, as in Figure [l] 
At low {a/R ^ 1) linear particle density, the optimal 
angular offsets that maximizes distance between neigh- 
boring particles is f2 = tt. However, when a/R drops 
below 2/\/3, the second neighbor in the axial coordi- 
nate (or "axial index") becomes the first neighbor in 
three-dimensional space, and so f2 = tt becomes unfa- 
vorable. As density increases further, the angles 27r/3 
and 47r/3 also become unfavorable due to third-neighbor 
interactions and so on, eventually generating a so-called 
Farey tree of unfavorable angles [7], where every new frac- 
tional multiple of 2-k is made by summing the numerators 
and denominators of the two previously adjacent multi- 
ples A growing plant evolves quasi-statically from 
one optimal to another as R/a increases, asymptoting 
to the golden angle f^i = 27r/ (r -\- 1) [r = (l \/5) /2] , 
ubiquitous in botany. Occasional "wrong turns" at bi- 




FIG. 1: A specimen of Mammillaria elongata displaying a 
helical morphology ubiquitous to nature, a magnetic cactus 
of dipoles on stacked bearings, and a schematic of a wrapped 
Bravais lattice showing the angular offset (screw angle) Q, and 
the axial separation a between particles. 

furcations during growth lead to alternative angles, the 
most common of which, flp — 2tt / (t +p) with p — 2,3, 
are called in botany second and third phyllotaxis. 

This appealing scenario requires experimental confir- 
mation, since it rests on the assumption that the ground 
states are helical Bravais lattices parameterized by the 
angular offset f2. Previous experimental reproductions 
of phyllotactic patterns did not interrogate the ground 
state of particles with long-ranged interactions; instead 
they examined kinetic processes in ferro-fluid droplets [5] 
or packings of hard spheres [9j . To properly test Levitov's 
model, we constructed two versions of a magnetic cac- 
tus consisting of 50 outward pointing, dipolar permanent 
magnets (spines) mounted on stacked coaxial bearings 
free to rotate about a vertical axis (stem) , as in Figure [T| 
The system can be annealed into a lower energy state by 
mechanical agitation [10]. We find that the final states so 
obtained consistently assume phyllotactic spirals whose 
angles precisely reproduce first {^i), second (ri2) and - 
for a cactus with longer spines/dipoles - third phyllotaxis 
(r^a), as expected in a biological cactus and as shown in 
Figure 2. Often, as also seen in Botany [11], the experi- 
mental system fragments into spirals of different fip, cor- 
responding to two or more local energy minima separated 
by a domain wall; this is expected given the degeneracy 
and low dimensionality. Numerical optimization via a 
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Axial Index 

FIG. 2: Measured angular ofltsets between successive mag- 
nets for magnetic cacti with short (top) and long (bottom) 
magnets, plotted versus axial index. Flat regions are perfect 
spirals and steps are kinks between domains of different spi- 
ralling angles. Dotted lines give the phyllotactic angles Qi, 
Q2, ^3 and 2n — ^1^ defined in the text. Dashed lines are 
minima of the magnetic lattice energy (insets) calculated by 
interpolating the measured pair- wise magnet-magnet interac- 
tion. The magnetic cactus recapitulates botany. 



structural genetic algorithm confirms and extends these 
results in axially unconstrained systems {lOj . 

Unlike its botanical analogue, the magnetic cactus 
can access dynamics. To better understand it, we 
first describe its fascinating self-similar energy land- 
scape [21 [7]. The essential results apply equally to any 
smooth moderate-ranged repulsive interaction. We con- 
sider dipoles p directed radially outward, interacting via 
Vt^j = Pi ■ Pj/rfj ~ 3{p^ ■ r,^j){pj ■ rjj)/rf The dipoles 
can rotate on their angular coordinate 9n but are con- 
strained axially to equal spacings z„ = an, n the axial 




FIG. 3: Top: Lattice energy V{n) (in units of e = p R~'^) 
versus screw angle for successively halving values ofa/R start- 
ing from 0.5, showing more highly degenerate curves. Bottom: 
energy curve for a/R = 0.15, with six distinct minima. The 
highest rank observed, J = 6, matches Eq. [l] 

ergy. To study dynamics, kinetic energy can be added: 
E = lIT,^^^'^ + Figure R plots the lattice en- 
ergy V{fl) versus angular offset IT for spiral lattices (i.e. 
On = ^ n) at various densities. As discussed above, 
the lattice energy develops nearly degenerate peaks cor- 
responding to commensurate spirals = 27r i/j with 
i,j coprime (i.e. having no common divisors): for these 
structures 6k — 0k+j , so is dominated by particles fac- 
ing each other at a distance ja; we call these peaks of 
rank j. For given a/R, there is a maximum rank J [10] 
(corresponding to the smallest peaks in Fig.[3|, 
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index. V = 1/2 X; 



is the total potential en- 



[[ ]] denoting the integer part. The minima between the 
peaks also become nearly degenerate with increasing den- 
sity, as each particle sees the others incommensurately 
averaged. Via number-theoretical considerations the de- 
generacy is found D = ^ J'^ + O {J log J) cx |f [TD]. 
A stable structure corresponding to a minimum brack- 
eted by peaks of rank ji and j2 is a spiraling lattice 
where nearest neighbors are at axial displacements ±aji , 
±0^2 and second nearest neighbors are at ±a(ji + J2) or 
±a(ji— j2) [3 . Thus ji andj2, which are coprime [7J, give 
the number of crossing helices (in Botany, parastichies) 
that cover the lattice by connecting nearest neighbors. 

This striking potential energy landscape returns non- 
monotonic dispersion relations in linearized dynamics, 
with many classical rotons and maxons. While rotons in 
superfluid "^He are associated with density fluctuations 
and generally held to be intrinsically quantum mechani- 
cal |12j . phyllotactic rotons have a new, purely geomet- 
rical and classical origin. Consider oscillations around 
a stable spiral Q.: On — ^ n + cos [kn — Lo{k)t\. The 
normal mode frequency is 

uj{kf = 2luI [1 - cos(jifc)] + 2u;l [1 - cos(j2fc)] + • • • (2) 
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FIG. 4: The phonon dispersion relations for the six phyl- 
lotactic lattices of Figure [Sj with multiple rotons and max- 
ons. Corresponding parastichy numbers {ji, j2) (defined in 
the text) and the three most strongly interacting neighbors 
[ill ia, is] (wiji > U2j2 > ijJsjs) are given. The simple es- 
timate for the number of rotons and maxons, 2ji — 1, holds 
for all but (2,5). Each spectrum is offset by 5r~^ for clarity, 
where r — e^^/^J^/^ is the unit of time [10) . 

where jVn is the axial index of the m-th neighbor and 
ujj^^ is the second derivative of wo,j„ about the equilib- 
rium spiral f2. For simplicity, let us truncate this ex- 
pansion at nearest neighbors. Crucially, nearest neigh- 
bors in the axial coordinate are not necessarily near- 
est neighbors in three dimensions, thus {ji, J27 js, ■■■}¥' 
{1,2,3,...}. As explained before, the nearest neighbors 
jij J2 are instead the parastichy numbers for fi, and thus 
coprime [3 : it follows that the dispersion relation of 
Eq. [2] truncated to just (ji,j2), has a periodicity in- 
terval [— TT, tt) (which is not a Brillouin zone, since our 
lattice is axially aperiodic), and is non-monotonic, but 
instead shows 2ji — 1 minima and maxima - or rotons 
and maxons (we order ji, j2 so that ujjj^ji > Wjjj^)- Fig- 
ure |4] shows the exact dispersion curves for each of the 
stable structures in the top panel of Figure [3] extrema 
follow the 2ji — 1 rule in all but one case. Physically, 
why are there rotons? Although increasing axial wave- 
vector drives particles nearby along the axis increasingly 
out of phase, it can drive particles that are neighbors in 
three dimensions more nearly in phase. This discrepancy 
between nearest neighbors in ID and 3D drives rotonic 
behavior. Simulations confirm the localization of energy 
and momentum that is expected for rotons lU]. 

While the linear dynamics of phyllotaxis provides ro- 
tons, the non/mear regime generates a new class of highly 
non-local topological solitons with unusually rich proper- 
ties. One-dimensional degenerate systems are generally 
entropically unstable against domain wall formation ^13] : 
we did see kinks both numerically and experimentally 
(Figure 2). Can these kinks travel as topological soli- 
tons? Experimentally, we observed the magnetic cactus 
expelling a higher energy domain by propagating its kink 
along the axis. Such axial motion of a kink between two 
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FIG. 5: Two phyllotactic solitons approach, collide and trans- 
form. Top: The screw angle plotted versus space and time 
for a collision of two phyllotactic solitons. Plateaux of dif- 
ferent height correspond to spirals of different angle. Middle: 
The potential energy V drops as the lower-energy domains 
advance, converting into rotational kinetic energy of the new- 
born domain (see online animation S5). Deviations starting 
around t = 2 arise from precursor elastic waves propagating 
in front of the solitons that ensure local angular momentum 
conservation [T3]. Bottom: The lattice energy V{Q) corre- 
sponds to the physically realized magnetic cactus of Figure 2 
(bottom), with minima at ni/27r = 0.23, Q.^jl-K = 0.28, 
Q,zl2-K — 0.38. For clarity, weak high-wave-vector phonons 
were smoothed by spatial averaging, here and in Figure [6] 

domains of different helical angles confronts a dilemma: 
helical phase is unwound from one domain at a different 
rate than it is wound up by the other. Numerical simula- 
tions show that the moving domain wall solves this prob- 
lem by placing adjacent domains into relative rotation 
(S5). If we compare this behaviour with the paradig- 
matic case of sine-Gordon-like 1-D topological solitons, 
which separate essentially equivalent static domains and 
can trave at any subsonic speed [13,, we see that phyl- 
lotactic domain walls instead separate regions of different 
dynamics: energy and angular momentum flow through 
the topological soliton as it moves, rather than being con- 
centrated in it, and its speed vk is tightly controlled by 
energy-momentum conservation, phase matching at the 
interface and boundary conditions. 

As a simple, symmetric example, consider a low- 
density system which reproduces our experimental cac- 
tus, prepared at rest in a metastable spiral at fli with 
free boundaries. As depicted in Figure |5] lower-energy 
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FIG. 6: The topological solitons can decay, merge, and change 
identity. An initial unstable domain boundary decays into 
two solitons. The left-hand soliton bounces off the free bound- 
ary, maintains its shape, propagates rightward faster than the 
other, collides with it, merges, and transforms into a third 
topological soliton with a different characteristic speed. S5 
provides an animation. The density is a/R — 0.15. 

domains of angle spontaneously form at the edges 
and move inward at a fixed speed while rotating uni- 
formly at 9 = vk (f^2 — ^i)- As the rigid f2i — domain 
wall advances, it converts the potential energy diS^erence 
V2 — Vi into rotational kinetic energy: / {fl2 — ^1)'^ v\t — 
2 (V2 — Vi) VKt, which requires a propagation speed of 
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The measured speed, 22.1 a s~^ (see Figure [5]), agrees 
well with the value predicted from Eq. [s] 23.4 as~^, 
the discrepancy likely coming from phonon radiation (ex- 
pected for a soliton on a discrete lattice) . 

An extraordinarily rich phenomenology, well beyond 
that seen in traditional soliton models [14j emerges from 
more complex numerical simulations: kinks of different 
species merge, decay, change identity upon collision, and 
decompose at high temperature into a sea of constituent 
lattice particles. Figure |6] depicts an unstable domain 
wall decaying into two topological solitons, one of which 
reflects off a free boundary then catches up and merges 
with the other to form a new soliton of different charac- 
teristic speed. A continuum analytical model that takes 
into account the rotational kinetic energy of the domains 
can predict and classify many of these phenomena |15j . 

Because dynamical phyllotaxis is purely geometrical in 
origin, this rich phenomenology of new excitations could 
appear across nearly every field of physics. Indeed phyl- 
lotactic domain walls have already been seen, but not 
recognized, in simulations [5 of cooled ion beams ^ 
where the system self-organizes into concentric cylindri- 
cal shells. Trapped ions or dipolar molecules could attain 
a J of several tens \A\ . Colloidal particles on a cylindri- 
cal substrate provide a highly damped version |16j . and 
polystyrene particles in air (as used to investigate |17j 
the KTHNY theory of 2D melting [18]) have reason- 
ably low damping and long-range interaction. With 



charged repulsive particles and the constraint to equal 
axial spacing released, each topological soliton carries a 
different charge ITS'. Strikingly, this charge is not con- 
served in soliton collisions, since it can be lent to or 
borrowed from a weakly pinned [19] lattice. This could 
happen in Wigner crystals on suspended, weakly doped 
semiconducting carbon nanotubes or boron nitride nan- 
otubes. Phyllotactic degeneracy requires small enough 
a/R = r1/2R^. Since > 37 glT, from Eq. [l]we have 
that R > f^Q-o — 5. 5 J A for a maximum peak rank J. 
flo is the Bohr radius for an appropriate effective mass 
and dielectric constant. If oq — 0.53 A, the minimal de- 
generacy, D — 2 and J = 2, occurs already for i? ~ 11 
A. Finally, the extension of this fascinating geometrical 
structure to the quantum regime is particularly intrigu- 
ing. 
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